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Description 

[REAL-TIME VELOCITY AND 
PORE-PRESSURE PREDICTION AHEAD 

OF DRILL BIT\ 

Background of Invention 

[0001] The invention relates generally to subsurface drilling op- 
erations and particularly to techniques for predicting pore 
pressures ahead of a drill bit while drilling a borehole 
through subsurface formations. 

[0002] a column of drilling fluid, usually referred to as "mud," is 
customarily provided in a borehole while drilling the bore- 
hole through subsurface formations. Usually, the weight 
of the mud is carefully selected such that the hydrostatic 
pressure gradient in at least the uncased section of the 
borehole is above the pore pressure gradient and below 
the fracture pressure gradient in the surrounding subsur- 
face formations. If hydrostatic pressure gradient is lower 
than pore pressure gradient, a kick or blowout may occur. 
If hydrostatic pressure gradient is higher than fracture 



pressure gradient, lost circulation may occur. Fracture 
pressure typically increases rapidly with depth so that 
maintaining hydrostatic pressure gradient below fracture 
pressure gradient after drilling an initial section of the 
borehole is usually less of a problem. Pore pressure on 
the other hand generally follows a less predictable pat- 
tern. To avoid drilling hazards, it is desirable to know the 
pore pressure gradient ahead of the drill bit so that the 
mud weight needed to provide the desired hydrostatic 
pressure gradient in the borehole can be determined prior 
to drilling an interval ahead of the drill bit. 
[0003] Methods are known in the art for predicting pore pres- 
sures ahead of a drill bit using well log data and/or seis- 
mic survey data. One common method of predicting pore 
pressures ahead of a drill bit using well log data involves 
determining a normal compaction curve from a well log, 
e.g., a sonic log, in combination with an appropriate geo- 
logical model. The normal compaction curve corresponds 
to the increase in formation density that would be ex- 
pected as a function of depth assuming absence of abnor- 
mal pressure. During drilling, a logging-while-drilling 
(LWD) sonic log is obtained and compared to the normal 
compaction curve. A consistent slowing trend of the LWD 



sonic log away from the normal compaction curve is used 
as a likely indicator of increased pore pressures ahead of 
the drill bit. This expected increase in pore pressures can 
be estimated from the amount of departure of the LWD 
sonic log from the normal compaction curve. The accuracy 
of this method largely depends on the accuracy of the well 
log data used in generating the normal compaction curve. 
Since this method relies on well logs, the well has to be 
drilled into the over-pressured zone for it to be detected. 
[0004] u.S. Patent 5,130,949 (issued to Kan et al.) discloses a 
method of predicting pore pressures ahead of a drill bit 
using well log data and surface seismic data. The method 
involves using well log data to estimate a shale fraction 
for a subsurface formation as a function of depth and the 
derivation of a shale compaction trend as a function of 
depth at intervals where the estimated shale fraction ex- 
ceeds a threshold. A translation curve expressing pore 
pressure gradient as a function of sonic interval transit 
time departure from the shale compaction trend is also 
derived using the well log data. Seismic observations are 
performed along a surface line, and interval transit times 
as a function of depth along the surface line are estimated 
from the seismic observations. The departures of the 



seismic interval transit times from the shale compaction 
trend as a function of depth for points along the surface 
line are computed. The departures are translated into 
pore pressure gradient predictions using the translation 
curve. The patent suggests improving pore pressure pre- 
dictions in deviated wells by adjusting interval transit time 
data for the shallower portion using check-shot data and 
applying the adjustments to an entire seismic section to 
obtain better depth. 
[0005] u.S. Patent 5,144,589 (issued to Hardage) discloses a 

method of estimating pore pressures ahead of a drill bit 
using drill-noise seismic. During drilling, the noise of the 
drill bit as it impacts the earth is used as a seismic source. 
Some seismic signals propagate directly from the drill bit 
to the surface. Some seismic signals propagate down- 
wardly and are reflected back to the surface. The direct 
signals are used to determine interval velocity for each 
formation through which the drill bit has drilled, and the 
interval velocity is continuously updated as the drill bit 
penetrates the earth. The reflected signals are used to de- 
termine acoustic impedance for each formation ahead of 
the drill bit. The interval velocity data and acoustic 
impedance information are combined to produce a log- 



like impedance estimation curve, which reflects the pore 
pressures ahead of the drill bit. In particular, the low fre- 
quency velocity trend immediately above the drill bit is 
extrapolated to produce the low frequency velocity behav- 
ior for a short distance, e.g., 100 to 500 feet, below the 
drill bit. This low frequency behavior is then used to cor- 
rect the acoustic impedance data for the formation ahead 
of the drill bit. Drill-noise seismic is generally inefficient 
in soft sediment, highly deviated boreholes, and while 
drilling with certain types of bits, such as polycrystalline 
diamond compact bits. Furthermore, estimation of acous- 
tic impedance ahead of the bit from reflection data rarely 
works and is not reliable. 
[0006] a new technique called seismic measurements- 
while-drilling (SMWD) has been applied to pore-pressure 
prediction ahead of a drill bit. In one implementation, a 
pore-pressure map is generated using surface seismic 
data. To obtain the map, velocities are estimated from 
pre-stack surface seismic data, and a velocity- 
to-pore-pressure transform appropriate for the area is 
used to convert the velocities to pore pressures. The 
pore-pressure map has the coordinates of the seismic 
shot position (horizontally) and seismic travel time 



(vertically). During drilling, SMWD is used to deliver real- 
time check-shot data. The check-shot data includes travel 
times of seismic waves generated at the surface as a func- 
tion of depth. The check-shot data are used to place the 
drill bit on the pore pressure map, thereby allowing the 
pore pressures ahead of the drill bit to be determined. 
The accuracy of this method depends on the accuracy of 
the pore-pressure map predicted from the surface seismic 
data. In general, the accuracy of the pore-pressure map 
predicted in this manner decreases as depth increases be- 
cause the velocities estimated from the surface seismic 
data become less accurate with increasing depth. 
[0007] Thus a need remains for more robust techniques for esti- 
mating pore pressures ahead of a drill bit. 
Summary of Invention 

[0008] The invention provides a method of estimating velocity 
ahead of a drill bit which comprises generating seismic 
waves at a surface from at least two different source posi- 
tions in the vicinity of a borehole, receiving seismic waves 
reflected from a reflector ahead of the drill bit at one or 
more locations in the borehole, determining travel times 
of the seismic waves received at the one or more locations 
in the borehole, and inverting the travel times to deter- 



mine a velocity of a formation ahead of the drill bit. In one 
embodiment, the method further comprises transforming 
the velocity into pore pressure of the formation. 

[0009] | n another aspect, the invention provides a system for es- 
timating velocity ahead of a drill bit which comprises at 
least one seismic source on a surface for generating seis- 
mic waves from at least two different source positions in 
the vicinity of a borehole and at least one seismic receiver 
in the borehole for detecting seismic waves reflected from 
a reflector ahead of the drill bit. The system further in- 
cludes a measurements-while-drilling telemetry system 
for transmitting data from the seismic receiver to the sur- 
face and at least one processor comprising instructions 
for determining velocities ahead of the drill bit using 
travel times of seismic waves reflected from the reflector. 
In one embodiment, the processor further comprises in- 
structions for converting velocities to pore pressures 
ahead of the drill bit. 

[0010] | n another aspect, the invention provides a method of es- 
timating pore pressure ahead of a drill bit which com- 
prises obtaining surface seismic survey data for a region 
of interest. During drilling of a borehole through the re- 
gion, the method further includes determining a travel 



time of a seismic wave generated from a surface of the re- 
gion to a location in the borehole when the drill bit is at 
selected depths in the borehole. The method further in- 
cludes determining a velocity from the travel time and the 
selected depths, inverting the surface seismic survey data 
to determine a velocity ahead of the drill bit while con- 
straining velocity between the surface and the drill bit to 
be consistent with the velocity determined from the travel 
time. In one embodiment, the method further comprises 
transforming the velocity ahead of the drill bit into pore 
pressure of a formation ahead of the drill bit. 
[001 1] other features and advantages of the invention will be ap- 
parent from the following description and the appended 
claims. 

Brief Description of Drawings 

[0012] Figure 1 is a flowchart illustrating a method of predicting 
pore pressures ahead of a drill bit according to an em- 
bodiment of the invention. 

[0013] Figure 2 illustrates a velocity model according to an em- 
bodiment of the invention. 

[0014] Figure 3A shows a drilling system suitable for practicing 
an embodiment of the invention. 

[0015] Figure 3B illustrates a SMWD system suitable for practicing 



an embodiment of the invention. 
[0016] Figure 4 is a flowchart illustrating a method of predicting 

pore pressures ahead of a drill bit according to another 

embodiment of the invention. 

[0017] Figure 5 illustrates a velocity model according to another 

embodiment of the invention. 
Detailed Description 

[0018] | n the following description, numerous specific details are 
set forth in order to provide a thorough understanding of 
the invention. It will be apparent, however, to one skilled 
in the art that the invention may be practiced without 
some or all of these specific details. In other instances, 
well-known features and/or process steps have been de- 
scribed in brief detail for clarity of the invention. 

[0019] Figure 1 is an overview of a method of predicting pore 

pressures ahead of a drill bit in accordance with one em- 
bodiment of the invention. The method can be used to 
predict pore pressure as the drill bit penetrates subsur- 
face formations. The method works as follows: when the 
drill bit is at a depth of interest in a borehole, one or more 
seismic sources on a surface are activated to generate 
seismic waves at multiple offset positions along the sur- 
face (100). The seismic waves are generated at predeter- 



mined times. In the borehole, one or more seismic re- 
ceivers on a drill string are used to detect seismic waves 
reflected off a target reflector ahead of the drill bit (102). 
Preferably there are multiple receivers or one receiver at 
multiple depths for velocity filtering of direct and reflected 
waves, as is well known in the art. The seismic waves re- 
ceived at the seismic receivers are processed, and the ar- 
rival times of the seismic waves are picked (104). The 
travel times of the reflected seismic waves are determined 
using the picked arrival times (106). The travel times of the 
reflected seismic waves are then inverted to determine ve- 
locities in formations ahead of the drill bit (108). The ve- 
locities ahead of the drill bit are transformed into pore 
pressures using known velocity-to-pore-pressure rela- 
tions (110), e.g., Hottman-Johnson relation. (Hottman, 
C.E., and Johnson, R.K., "Estimation of Formation Pres- 
sures from Log-Derived Shale Properties," JPT, Vol. 17, 
June 1965, pp 717-22.) 
[0020] | n s tep ioo, the seismic waves could be generated by a 

single seismic source that is offset from the borehole in a 
predetermined manner. Alternately, multiple seismic 
sources could be used to the generate the seismic waves 
at the multiple offset positions. The offset positions at 



which the seismic waves are generated relative to the 
borehole can range from zero to a maximum value. Typi- 
cally, the maximum offset is not greater than half of the 
depth of the drill bit. In general, the maximum offset de- 
pends on the depth of the drill bit, the depth of the target 
reflector, and general velocity trends. The number of off- 
set positions needed depends on the assumed velocity 
model of the zone between the drill bit and the reflector. 
For a simplistic two-velocity model, two offset positions, 
e.g., one at the borehole and the other at the maximum 
offset, are the minimum required. In practice, more offset 
positions may be required to combat noise and uncertain- 
ties in anisotropy and to account for complicated velocity 
models. 

[0021] Figure 2 illustrates a simple two-velocity model where 

subsurface formations 200 are viewed as consisting of two 
layers 202, 204. The layer 202 extends from the surface 206 
to the depth of a drill bit 208, which is located in a bore- 
hole 210 in the subsurface formations 200. The layer 204 
extends from the depth of the drill bit 208 to a target re- 
flector 212 ahead of the drill bit 208. The velocities in the 
layers 202, 204 can be represented with a single average 
velocity respectively. For discussion purposes, let Vo be 



the average velocity from the surface 206 to the depth of 
the drill bit 208 and Vp be the average velocity from the 
drill bit 208 to the reflector 212. Further, let H be the offset 
position relative to the borehole 210 at which a seismic 
wave is generated, B be the depth of the drill bit 208, and 
R be the distance from the drill bit 208 to the reflector 212. 
If Snell's law is ignored (i.e., using straight ray paths), then 
the reflected-wave travel time, i.e., the time it takes a 
seismic wave generated from a source 214 at the surface 
206 to travel down to the reflector 212 and back up to a 
downhole seismic receiver 216, can be expressed as: 



^ F * J Vp Vo ' 2R + B w 



[0022] h, B, and Vo are known, and R and Vp are unknown. Vo 
can be obtained from a combination of pre-drill informa- 
tion, well-log data, LWD logs, and real-time check-shot 
data. Because there are two unknowns, at least two reflec- 
tion events, producing two reflected-wave travel times, 
are needed to determine R and Vp. Two reflection events 
can be obtained by generating two seismic waves 218, 220 
at two offset positions at the surface 206 and measuring 
the travel time of the ensuing reflected waves 222, 224, re- 



spectively. The solutions would most likely be found using 
numerical techniques. Therefore, the more the travel time 
data available, the better the results can be. 
[0023] a more complicated velocity model can be used to im- 
prove the accuracy of the method described above. For 
example, the subsurface formations 200 can be treated as 
consisting of several horizontal layers instead of just two 
layers. In this case, a series of velocities would be defined 
for the section between the surface 206 and depth of the 
drill bit 208. The velocities of the horizontal layers in this 
section could be obtained from real-time check-shot data 
and/or well log data, e.g., sonic logs. The real-time 
check-shot data is obtained by generating seismic waves 
at the surface. The seismic waves are generated in the 
vicinity of the borehole for selected bit depths, usually 
corresponding to depths at which drill pipe connections 
are made. The seismic waves (direct arrivals) are detected 
downhole by one or more seismic receivers on a drill 
string. The detected waveforms or the first arrival times of 
the detected waveforms are then sent to the surface via 
measurements-while-drilling (MWD) telemetry when mud 
circulation resumes. Real-time check-shot data collected 
for two different bit depths can be used to compute the 



interval velocity between the two bit depths. 

[0024] The section between the drill bit 208 and the target reflec- 
tor 212 may be modeled as one effective layer having an 
average velocity or may be represented by a parameter- 
ized curve that allows for some smooth variation in veloc- 
ity with depth. In this case, the reflected -wave travel time 
could be computed by a ray tracing algorithm, and the 
unknown parameters, i.e., distance of drill bit to reflector 
and velocities ahead of drill bit, can be estimated by fit- 
ting the computed reflected -wave travel times to the 
measured reflected -wave travel times, for example, using 
least-squares technique. 

[0025] Figure 3A shows an example of a drilling system 300 suit- 
able for practicing the present invention. The drilling sys- 
tem 300 includes a rig 302 positioned on a surface 304. 
The rig 302 is shown as a land rig, but could also be an 
offshore rig in other embodiments. The rig 302 includes a 
derrick 305 with an attached drill string 306. The drill 
string 306 extends into a borehole 308 through a wellhead 
310. A drill bit 312 is appended to the drill string 306 for 
drilling the borehole 308 through subsurface formations 
314. As the drill bit 312 penetrates the subsurface forma- 
tions 314, mud pumped down the drill string 306 is forced 



out of the nozzles of the drill bit 312 into the bottom of 
the borehole 308. The mud at the bottom of the borehole 
308 rises up the borehole 308 and is diverted to a mud re- 
turn system (not shown) on the rig 302. The drill bit 312 is 
located above a zone 318 in the subsurface formations 
314. The zone 318 may have abnormal pore pressures, for 
example. Therefore, it is helpful to know the pore pres- 
sures likely to be encountered in the zone 318 prior to 
drilling through the zone 318. This information is useful in 
selecting the weight of the mud needed to maintain the 
proper hydrostatic pressure gradient in the borehole 308 
while drilling the zone 318. 
[0026] | n one embodiment, the zone 318 includes a strong reflec- 
tor 320 or is located between the drill bit 312 and a strong 
reflector 320. The drill string 306 includes a downhole tool 
322, which includes one or more seismic receivers 324 that 
can detect seismic waves reflected off the reflector 320 as 
well as other seismic waves propagating in the subsurface 
formations 314. The seismic waves reflected off the reflec- 
tor 320 are generated by one or more seismic sources 326 
on the surface 304. The seismic sources 326 are offset 
from the borehole 308, where the offset value can range 
from zero to a maximum value depending on the depth of 



the drill bit 312, the depth of the reflector 320, and general 
velocity trends. The number of seismic sources 326 or off- 
set positions would depend on the assumed velocity 
model of the zone between the drill bit 312 and the reflec- 
tor 320. It should be noted that one or more seismic 
sources may be moved between the various offset posi- 
tions, i.e., it is not necessary to have a seismic source per 
offset position. 

[0027] Also, it is not necessary to have the drill bit 312 in the 
same depth position for various offset positions of the 
source. The drill bit 312, and the seismic receivers 324, can 
be anywhere above the look-ahead zone 318 for different 
offset positions as long as they can detect the reflected 
events from horizon 320. In this case, numerical ray trac- 
ing algorithms can be used to compute travel times of re- 
flected events. The unknown velocities ahead of the bit 
can be estimated by fitting the computed reflected wave 
travel times to the measured reflected wave travel times. 
For example, this can be accomplished using least 
squares techniques. This may be important in practice 
since one or several sources can be moved between offset 
positions as the drilling proceeds above zone 318. 

[0028] Figure 3B illustrates a SMWD system that includes the 



seismic receivers 324 and the seismic sources 326. The 
seismic sources 326 are coupled to a firing system 328, 
which is coupled to a processor 330 and clock 332. The fir- 
ing system 328 activates the seismic sources 326 in re- 
sponse to instructions from the processor 330. The pro- 
cessor 330 is programmed to cause the firing system 328 
to activate the seismic sources 326 according to a sched- 
ule known by the seismic receivers 324. A telemetry sys- 
tem 334 is also coupled to the processor 330 and receives 
data measured by the seismic receivers 324. The seismic 
receivers 324 are coupled to a signal processor 336. Mem- 
ory 338 and clock 340 are also coupled to the signal pro- 
cessor 336 so that recordings may be made of waves de- 
tected by the seismic receivers 324 in synchronization with 
firing of the seismic sources 326. A motion sensor 342, a 
mud flow sensor 344, and MWD telemetry equipment 346 
are coupled to the signal processor 336 for transmitting 
data to the telemetry system 334. 
[0029] Returning to Figure 3A, when it is desired to predict pore 
pressures ahead of the drill bit 312, circulation of mud and 
motion of drill string 306 are preferably stopped. There- 
after, the seismic sources 326 are activated in a predeter- 
mined sequence known to the downhole tool 322 to gen- 



erate seismic waves. The seismic waves propagate down 
the subsurface formations 314. The seismic waves reach- 
ing the reflector 320 are reflected back towards the sur- 
face 304. The seismic receivers 324 detect the up-going 
waves, i.e., the reflected waves, as well as the down-going 
waves, i.e., the direct waves. The detected waveforms (or 
representations of the detected waves) may be stored in 
memory (338 in Figure 3B). 
[0030] | n one embodiment, when mud circulation resumes, the 
detected waveforms stored in memory (338 in Figure 3B) 
are sent to the surface via MWD telemetry (346 in Figure 
3B) and received at the surface via telemetry system (334 
in Figure 3B). The detected waveforms may be sent to the 
surface in full or compressed format. The time required to 
send the detected waveforms depends on the available 
MWD telemetry rate as well as the number of bits/sample 
required, which is related to the strength of the detected 
waves. In this embodiment, the processor (330 in Figure 
3B) at the surface runs algorithms to separate the up- 
going (reflected) waves from the down-going (direct) 
waves. Techniques for making this separation are known 
in the art. The reflected waves are then processed to pick 
their arrival times. Techniques for processing waveforms 



to pick their arrival times are known in the art. 
[0031] | n another embodiment, instead of sending the detected 
waveforms to the surface for processing, the processor ( 
336 in Figure 3B) may include instructions for separating 
the up-going (reflected) waves from the down-going 
(direct) waves. The reflected waves may be processed 
downhole for picking of their arrival times or the reflected 
waves may be sent to the surface for processing. When 
mud circulation resumes, the arrival times of the reflected 
waves or the reflected waves themselves are then trans- 
mitted to the surface via MWD telemetry (346 in Figure 3B) 
and received at the surface via telemetry system (334 in 
Figure 3B). In the latter case, further processing at the 
surface is needed to pick the arrival times of the reflected 
waves. 

[0032] Once the arrival times of the reflected waves are known, 
the travel times of the reflected waves can be determined 
and used to compute the velocities ahead of the drill bit. 
The velocities can then be transformed into pore pres- 
sures ahead of the drill bit. In operation, the clock (332 in 
Figure 3B) at the surface keeps track of the time the seis- 
mic waves are generated, while the clock (340 in Figure 3B) 
downhole keeps track of the time the reflected waves are 



detected. These two clocks allow the arrival times and 
travel times of the seismic waves to be determined and al- 
low coordination of the generation and detection of the 
seismic waves. It is preferable that both clocks are in sync 
and are accurate enough to allow velocities to be com- 
puted with desired accuracy. The desired accuracy of the 
clocks depends on the particular method used in comput- 
ing the velocities. 
[0033] Figure 4 is an overview of a method of predicting pore 
pressures ahead of a drill bit according to another em- 
bodiment of the invention. The method starts by acquiring 
surface seismic data for the geological area of interest ( 
400). Methods for acquiring surface seismic data are well 
known in the art. Briefly, surface seismic data can be ac- 
quired by shooting multiple seismic waves at different lo- 
cations on the surface of the geological area. As an exam- 
ple, the seismic sources could be explosive charges (on 
land) or air guns (at sea). The seismic waves travel down 
the subsurface formations and pass through or reflect off 
a target. In reflection seismic surveys, the reflected seis- 
mic waves are detected and recorded using, for example, 
geophone (on land) or hydrophones (at sea). The seismic 
waves are recorded at the surface at various distances 



from the seismic shot positions. For each source-receiver 
pair, the travel time of the seismic wave from the source 
to the receiver is determined from the arrival time of the 
recorded seismic wave. 
[0034] As well known in the art, the surface seismic data can be 
used to estimate velocities of the subsurface formations 
(see, for example, Ozdogan Yilmaz, "Seismic Data Pro- 
cessing, " Society of Exploration Geophysics, 1987). These 
estimated velocities are typically known as pre-drill veloc- 
ities because the surface seismic data are usually obtained 
prior to actually drilling the borehole. A pore-pressure 
map can be generated using these pre-drill velocities, al- 
lowing pore pressures ahead of the drill bit to be deter- 
mined as the borehole is drilled. However, the accuracy of 
the pre-drill velocities decrease with increasing depth, 
leading to less accurate pore-pressure prediction with in- 
creasing depth. To improve accuracy of pore-pressure 
predictions made from pre-drill velocities, the inventors 
propose updating the surface seismic data with real-time 
check-shot survey data. With real-time check-shot survey 
data, accurate measurements of velocities from the sur- 
face to specific depths in the borehole, between the sur- 
face and the drill bit, can be obtained. These accurate ve- 



locities can then be used in place of the pre-drill velocities 
in the section from the surface to the drill bit while pre- 
dicting pore pressures ahead of the drill bit. 
[0035] From the foregoing, the method further includes conduct- 
ing check-shot surveys as the borehole is drilled and de- 
livering the check-shot information to the surface in real- 
time (402). A check-shot survey involves firing a seismic 
source at the surface, near the opening of the borehole. 
The seismic wave generated travels down the subsurface 
formations and is detected and recorded by a receiver in 
the borehole. This is done when the receiver, which is 
typically placed on the drill string, is located at the depth 
of interest. The drill string may include a downhole tool 
that can process the recorded wave to determine the ar- 
rival time of the wave at the receiver location. Alterna- 
tively, processing can be done at the surface. Either the 
recorded waveform or the picked arrival time is sent to 
the surface via MWD telemetry. To reduce noise in the 
data, it is preferable to temporarily stop drilling while 
making the check-shot survey. The waveform or arrival 
time can then be transmitted to the surface when drilling 
resumes. To avoid lost drilling time, it is preferable to 
make check-shot measurements during pipe connections. 



The travel time of the seismic wave from the seismic 
source to the receiver can be determined from the arrival 
time. The travel time can then be used to calculate aver- 
age velocity of the subsurface formations above the drill 
bit. Sometimes, it is convenient and accurate to subdivide 
the subsurface formations into a series of intervals. In this 
case, a series of travel times acquired in the manner de- 
scribed above can be used to calculate the average veloc- 
ity for each interval. 
[0036] The next step is to re-invert the surface seismic data us- 
ing the check-shot data to obtain improved velocities 
ahead of the drill bit (404). Figure 5 shows a simple veloc- 
ity model that can be used in re-inverting the surface 
seismic data. In this model, the drill bit 500 is shown at a 
desired position in the borehole 502, above a zone of in- 
terest 518, which is located above or includes a target re- 
flector 506. Multiple seismic shot positions 508 and re- 
ceiver positions 510 are superimposed on the model. The 
seismic shot positions and receiver positions are used 
prior to drilling the borehole to acquire surface seismic 
data, as previously explained. In this case, the seismic 
shot positions and receiver positions of interest are those 
that measure reflection events in the vicinity of the bore- 



hole 502 or have a common midpoint at the borehole 502. 
Also shown in the figure is the seismic shot position 512 
used in making check-shot survey. The model assumes 
that the subsurface formations 514 can be viewed as a 
two-layer system, where one layer 516 is defined above 
the depth of the drill bit 500 and the other layer 518 is de- 
fined between the depth of the drill bit 500 and the reflec- 
tor 506. 

[0037] Assuming that Snell's law can be ignored (i.e., using 

straight ray paths), then using the model described above, 
the travel time observed by a source-receiver pair on the 
surface can be expressed as: 



y J Vp Vo R + B w 



[0038] where H is the offset distance between the source and the 
receiver, B is the bit depth, Vo is the average velocity from 
the surface to the bit depth, R is the distance from the 
drill bit to the reflector, and Vp is the average velocity in 
the zone between the drill bit and the reflector. H and B 
are known. T can be obtained for various values from the 
surface seismic data. Returning to Figure 4, Vo is com- 
puted from the real-time check-shot survey in step (402) 



instead of from the surface seismic data obtained in step ( 
400). R and Vp are unknown. To obtain R and Vp, equation 
(2) above is inverted. At least two sets of T and H values 
for the selected bit depth B are needed to compute R and 
Vp, but accuracy generally improves with more data sets. 
Vo is the same for a given bit depth B for all source offset 
positions. 

[0039] Once the velocity ahead of the drill bit (Vp) is known, the 
pore pressure ahead of the drill bit can be determined us- 
ing any suitable velocity-to-pore-pressure relation, e.g., 
Hottman-Johnson relation, (406). To improve the accuracy 
of the pore-prediction method, a velocity model that is 
more complicated than the one described above can be 
used. For example, instead of treating the subsurface for- 
mations as a two-layer system, the subsurface formations 
could be divided into several horizontal layers, each layer 
having its own associated average or interval velocity. The 
interval velocities for the layers between the surface and 
the drill bit can be obtained from the real-time check- 
shot data. This allows variations in the velocities above 
the drill bit to be accounted for while determining the ve- 
locities ahead of the drill bit. The section between the drill 
bit and the reflector may be modeled as one effective 



layer having an average velocity or may be represented by 
a parameterized curve that allows for some variation with 
depth. 

[0040] it w i|| be apparent to those skilled in the art that this in- 
vention may be implemented by programming one or 
more suitable general-purpose computers having appro- 
priate hardware. The programming may be accomplished 
through the use of one or more program storage devices 
readable by the computer processor and encoding one or 
more programs of instructions executable by the com- 
puter for performing the operations described above. The 
program storage device may take the form of, e.g., one or 
more floppy disks; a CD ROM or other optical disk; a 
magnetic tape; a read-only memory chip (ROM); and other 
forms of the kind well known in the art or subsequently 
developed. The program of instructions may be "object 
code," i.e., in binary form that is executable more-or-less 
directly by the computer; in "source code" that requires 
compilation or interpretation before execution; or in some 
intermediate form such as partially compiled code. The 
precise forms of the program storage device and of the 
encoding of instructions are immaterial here. Thus these 
processing means may be implemented in the surface 



equipment, in the tool, or shared by the two as known in 
the art. 

[0041] The invention provides several advantages over conven- 
tional techniques. The method of predicting pore pres- 
sures ahead of a drill bit is based on getting look-ahead 
velocities in near real-time and relating these to look- 
ahead pore pressures by known transforms. The method 
is superior to techniques using only surface seismic ve- 
locities because it incorporates accurate real-time mea- 
surements of velocities measured above the drill bit to in- 
crease the accuracy of estimated velocities ahead of the 
drill bit. Velocities from surface seismic alone degrade, 
i.e., become less accurate as depth increases. Further, 
methods that invert impedance profiles ahead of the drill 
bit are known to have limited success. One reason is that 
reflection data is a good measure of sharp impedance 
variations ahead of the drill but, but it does not contain 
"slow" impedance variations. Reflection data does not give 
the "velocity" itself ahead of the bit, instead it gives 
"change in velocity" or the derivative of the velocity. 
Therefore, one needs to "assume" the velocity trends 
ahead of the bit to use these techniques. The present in- 
vention, on the other hand, uses the actual measured ve- 



locities ahead of the bit. 
[0042] while the invention has been described with respect to a 
limited number of embodiments, those skilled in the art, 
having benefit of this disclosure, will appreciate that other 
embodiments can be devised which do not depart from 
the scope of the invention as disclosed herein. For the 
purposes of this specification it will be clearly understood 
that the word "comprising" means "including but not lim- 
ited to", and that the word "comprises" has a correspond- 
ing meaning. 



